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ABSTRACT 


\ 

\ 



Electrostatic plugging of multidipole cusps is investigated. 

The plugging is produced by placing positively biased electrodes on 
both sides of all magnetic cusps. Experimental observations show 
the plasma potential follows the electrode potential and the plasma 
electron density increases by a factor of about two as the electrode 
bias increases from zero. The addition of a negatively biased grid 
at one end of the chamber does not change the response of the plasma 
electron density to the electrode bias. A model, which considers 
plasma losses along the magnetic cusps as well as by diffusion trans¬ 
verse to the field lines, is shown to qualitatively agree with the 
data. The elevated interior plasma potential permits the identifi¬ 
cation and investigation of secondary electron production at the 
walls. Data is presented which shows that secondary electron emission 
rate at the chamber walls can be comparable to plasma production rate. 
A measurement of the relative ion densities oif Hydrogen is attempted, 
but the results are inconclusive as to the effect improved confine¬ 
ment has on the various ion species densities. 
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I. INTRODUCTION 


The use of fields to confine laboratory plasma has been 
studied for many years. The use of magnetic fields, 1 to increase 
particle path lengths, is limited by losses along the field lines. 

For example, the use of cusps and surface magnetic fields has shown 
that the loss area for ions is bounded by UrgL, (Ref. 3) and 2r^L c 
(Ref. 4) where r H = , r^ is the ion gyroradius, r g is the 

electron gyroradius, and L q is the length of the magnetic cusps. 

2 

Attempts using electrostatic fields, to keep either ions or electrons 
from the chamber walls, have failed because of heating of the grid 
wires or capacitor plates, used to produce the fields, from plasma 
bombardment and also because of arcing problems. 

To improve plasma confinement, various authors have tried a 
combination of electrostatic and magnetic confinement. Dolan^ 
describes this technique as "... magnetic shielding of the grid of an 
electrostatic plasma confinement device or electrostatic plugging of 
a magnetic confinement device." Most of this early work has involved 
ring spindle cusp devices, with electrodes biased negative (with 
respect to the plasma). The present work is the first attempt at 
using electrostatic confinement in conjunction with the full line 
cusp surface magnetic fields. 
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This work vises a standard kO liter multi-dipole soup pot in a 
full line cusp configuration,^" and uses magnetically shielded 
electrodes to alter the internal plasma potential structure in a 
way that improves the plasma confinement. 

The apparatus used for this work is described in Section II. 
This section also describes briefly some of the diagnostic devices 
used, further descriptions being provided as needed. Section II 
also discusses the basic operating regime of the device. 

Experimental results are presented in Section III. Included 
are variations of plasma potential, density, and temperature with 
electrode bias. Also presented are results obtained with the device 
configured to simulate an ion source. 

A simple theory is presented in Section IV which gives good 
qualitative agreement with the results of Section III. 

Section V deals with secondary electron emission from the 
chamber wall due to electron bombardment. The appearance of the 
secondary electrons is a result of the way the plasma responds to 
the biased electrodes when the neutral pressure is less than 
10 ^ Torr. 

Section VI is devoted to seme experiments using hydrogen as 
the neutral gas. Hydrogen was used to determine the effect the 
electrodes had on the ion composition. 

Section VTI is a final summary of the work presented. 







II. APPARATUS 


The chamber used for this work was a 40 liter multi-dipole 
soup pot, with the permanent magnets in a full line cusp arrange - 
ment. The magnets produced a maximum field of 1 kilogauss at the 
poles. 

The basic device was modified by placing positively biased 
(with respect to the chamber wall) electrodes straddling each row of 
magnets^ (Fig. l). The electrodes were approximately 0.7 cm above 
the wall. At this distance the magnetic field was of order 500 Gauss, 
so the electrodes were magnetically shielded from the plasma. It 
was possible to change the gap width from 0.4 cm to 1.3 cm. This 
separation is on the order of the ion-gyroradius and much greater 
than the plasma electron gyroradius (r^ s 1 cm for argon, helium, 
and hydrogen; r g s 5 X 10 ^ cm). 

The plasma was produced by electron emission from negatively 
biased filaments located in the field free interior. The energy of 
these ionizing electrons (hereafter called primaries) was equal to 
the difference of the filament bias potential and the plasma 
potential <p p . The device was operated with low electron injection 
(discharge) currents (i^ £ 1A) which produces densities on the order 
of 10 10 cm”^, and low fractional ionization ( s 1%). 






The primary diagnostic tools used were collecting and emitt¬ 
ing Langmuir probes. The collecting probes, l/4" diameter Tantulum 
discs, were used to measure plasma electron temperature and density, 
and <Pp in the field free regions. The emissive probe (3 mm long, 
0.003 cm diameter, hot tungsten wire) was used to map the spatial 
plasma potential variation in the device. How we used the probes 
will be explained more in Section III and Section V. 

This work also required the use of mass spectrometers mounted 
inside and outside the chamber, and a gridded disc for launching 
pseudo-waves. This apparatus will be discussed more fully in 
Section VI. 


Chamber Operation 

The primaries are confined by the surface magnetic field. 
Previous experiments have shown that the leak width of the primaries 
is only the order of their gyroradius (0.03 cm for the maximum |b| 
of 1 kilogauss)Plasma electrons and ions are confined by a 
combination of magn etic fields and self-consistent electric fields, 
which for large electrode gaps (G > r^) are essentially the same as 
those in multidipole devices without electrodes. 

Consider operation with the electrodes biased at cp E such that 
ep E /T e » 1. For large gaps the plasma potential will be on the 
order of T g /e above ground, which is the anode. As the gap narrows 
to less than the self-consistent ion leak width the ions will be 


electrostatically confined at the electrodes and the ion leak width 








will be reduced to the gap width or slightly narrower. At the same 
time the electron leak rate will increase because of the presence of 
the positively biased electrodes. Improved ion confinement and 
reduced electron confinement will lead to an increase in the plasma 
potential until the electron and ion loss rates balance. As the 
plasma potential increases, electron losses to the wall are greatly 
reduced (exponentially as e eCp p/ T e) and electron losses to the posi¬ 
tive electrodes are increased. Thus the electrodes assume the 
function of the anode in the plasma and the plasma potential self- 
consistently adjusts to be within a few T g /e of the electrode 
potential. 

When * <p£ » T Q /e, plasma electrons are electrostatically 
confined at the cusps, as mentioned above, and plasma electron leaks 
are then dominated by diffusion across the between cusp magnetic 
fields to the electrodes. Primaries can still be lost at the cusps 
and this loss flux has an additional consequence. Cold secondary 
electrons with energies the order of a few eV are produced by the 
primaries at the walls and secondary electron currents can be com¬ 
parable to the current from the primaries.^ In ordinary multidipole 
devices these secondary electrons are indistinguishable from the 
plasma electrons. In this chamber they are distinguishable because 
the secondaries are accelerated through » 1 before reaching 

the interior of the device. These essentially monoenergetic electrons 
are readily confined by the surface magnetic fields, become spatially 










isotropic and are readily identified as an additional electron 

species. See Section V for more details. 

The energy of the primaries E is equal to e(cp - ) because 

P P F 

the plasma serves as a virtual anode for the filaments. Since the 

plasma potential <P^ generally follows cp^ (and q> F < 0), cp^ was varied 

to keep V - «p_ = constant. The ionizing electron energy E remains 
P * P 

essentially constant as tp £ is varied. Although the ionization cross 
section peaks below 100 eV (e.g., 90 eV for ergon); 10 inelastic 
collisions, multiple ionizations by single electrons, and other 
factors combine to give a maximum plasma density for E^ greater than 
100 eV, and this optimum E^ increase with increasing neutral 
pressure. 11 Above this energy the density stays close to the maximum 
value. For this work the primary electron energy E^ was chosen 
somewhat greater than the value for maximum plasma density so any 
residual small change in E^ does not have the effect of changing 
the ionization efficiency. 







Ill. EXPERIMENTAL RESULTS 


This section will present the various experimentally deter¬ 
mined. plasma parameters, plasma potential (cp^), electron temperature 
(T g ), and density (n) and their variation with electrode bias (cp £ ), 
neutral pressure (p Q ), and electrode separation (G). The operation 
was confined to two separate modes: l) multidipole plasma chamber, 
2 ) simulated ion source. 

The diagnostics used were emitting and collecting Langmuir 
probes. The collecting probe was typically placed near the center 
of the chamber. A voltage was swept across the collecting probe, 
and an I vs. V Langmuir probe characteristic was generated and 
recorded on an x - y plotter (Fig. 2). The method of getting 
parameters from these characteristics involved removing the primary 
electron portion by drawing the straight line through that part of 
the trace (Fig. 3). That line is the baseline, the zero current 
line for plasma electrons (secondary electrons are only observed at 
pressures lower than those used in this part of the experiment and 
will be ignored for present), and the natural logarithm of the 
plasma electron current is plotted against the applied voltage 
(Fig. 4). This semi-log plot will ideally have two parts. First, 
a linear part up to tp^, and a second part that is not collinear with 






the first. The voltage •where this discontinuity occurs is called 
the plasma potential. The voltage difference over which I changes 
by a factor e is the electron temperature T . 


There is a third, parameter measured on the characteristic and 
that is the plasma saturation current, the value of I at the dis¬ 
continuity mentioned above. This current can be written^ 


1/2 


= e n A (- 


( 1 ) 


A is the probe area, and T g is measured in energy units. For a 
1/4" diameter probe the density can be written 


n - 5.8 X 10 8 I^/JT (2) 

-3 * 

where n is in cm , is in mA, and T g is in eV. Now we have 
measured, or calculated the plasma parameters from the character¬ 
istic. 


Multi-Pinole Operation 

First consider the variation of <P p with applied electrode 
voltage «P E (Fig. 5). Note that generally follows <P £ , although 
initially > «P E ; then after ecp^/Tg > 2 the opposite is true. The 
difference between the applied voltage and plasma potential is on the 
order of T g . 
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T g is presented (Fig. 6) as a function of cp^. T g is shown to 
he a fairly constant function of for a specific parameter set, 
the error bars on T are of the order l/2 eV. Later it will be 
assumed (Section TV) that T g is constant and these data show that 
this assumption is warranted. 

The density also increases as cp £ increases. We call the ratio 
of the density to its value with zero electrode bias the enhancement 
ratio (7). Fig. (7) shows the variation of 7 as increases. 
Presented also in this figure is the dependence on p Q and gap width. 
Note the largest ratio of y is obtained when G = 12 mm and 

-<5 _k 

p = 2 X 10 . The lowest value of y is obtained for p = 1 x 10 

0 o 

and G = 6 mm. These results for y can qualitatively be compared to 
calculated values of 7. This is done in Section IV. 

Recall this experiment was initially intended to prove that 
the electrodes improve the ion confinement, but ions have not been 
discussed yet. We infer increased ion confinement from the electron 
density increase. Physically, what is impeding the ion loss? If 
the ion loss area is reduced then the electrodes are improving the 
ion confinement. 

Looking for the reduction in loss area required making equi- 
plasma potential contours around the cusp area. Since this area was 

magnetic, a collecting Langmuir probe couldn't be used. However, 

12 


Smith et al, 


showed that an emissive La n g m uir probe could give the 
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desired information. We used the inflection point method they 
present. 

The results are presented in Fig. (8). These results do 

show the formation of an ion trap. In this case the ion trap 

reaches 8 volts, but since the ions are very cold a trap of order 1 

volt would be sufficient in these devices. 

This figure also shows the reduction of the ion loss area. 

The loss area is defined as A_ = W L where L is the length of the 
L L c c 

magnetic cusp. W^ is an empirical loss width through which ions 

with energy e(<P ) + T, can flow out. («p ) means the plasma 

^ c *c 

potential in the center of the device. The upper figure shows this 

energy is s 0.5 + ^ and the lower figure shows Eg = 42 + T^ 

The corresponding loss widths are W = 9 mm and W = 5 mm. This 
L 1 L 2 

analysis was performed for many <P E , and a linear relationship was 
found between n and (W^)’ 1 (Fig. 9). 

Simulated Ion Source Operation 
Operating the chamber as a simulated ion source was done by 
covering one end of the cylinder by a grid (or a plate) and removing 
the magnets and electrodes. We again measured cp^, n, and T g as 
functions of <p £ , p Q . Gap width was set at 12 mm. 

Representative data are shown in Fig. (10). The points to 
notice are: l) <p^ is slightly less than «p£ at high electrode bias 
and insensitive to whether the grid is floated or grounded, 2) the 
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density does saturate, but the saturation occurs at higher cp r than 
when the grid is not present, and 3) the factor by which the density 
increases is greater for the grounded grid case. A simple calcula¬ 
tion explaining the variation of density for a grounded and floating 
grid is presented in the next section. 

T g , for the case shown in Fig. (lO), is graphed in Fig. (ll). 
Again, T g is a fairly constant function of «P E * 

Thus, overall the grid has very little impact on the varia¬ 
tion of the measured parameters of the plasma with electrode bias. 

Of course, the density is reduced because the grid presents a large 
loss area for ions and electrons, there being no magnetic shielding 
of the grid. However, the main point is that whether the chamber is 
operated as an ion-source or multidipole soup pot, the results are 
qualitatively the same (compare Fig. 7 and Fig. 10). 








IV. THEORY 


This section will present a simple model that predicts a 
density increase in qualitative agreement with the experiment. A 
basic assumption used at this time is that the plasma potential is 
increased as the electrode bias is increased, and about equal. In 
the second part of section IV this assumption will be proved. 

The quantities to consider are the electron loss rate along 
the cusp at low electrode bias, and the electron loss rate across 
the magnetic field at high electrode bias. These quantities are 
considered because they are the dominant electron loss mechanisms. 

The electron loss rate out the cusp can be expressed as: 


n v - e cp /T 

: -S-£ w L e p e 


(3) 


The ion loss rate can be written: 


■ -T 


GO 


j + /e is a Bohm loss flux which equals nc g /2. is the ion loss 

area at the cusps and equals p U P is an effective ion leak which 
has been shown^ for similar devices to equal Wr g r^, L c is the 
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total cusp length. A f ^, the area of the filaments, is only 0.1$ 
of Ayr, so we will ignore it. n g is the total primary and plasma 
electron density, but the primary density is small so we approximate 
n g by n, the interior plasma density. v g is the electron thermal 
velocity (v g = v'(T e /M g ), and c g is the ion-acoustic velocity 
(c g = . The exponential takes into account electron electro¬ 

static confinement along the cusp. The only quantity undefined is 

w which is an effective electron loss width, and the unknown of the 
c 

present calculation. 

Equating the electron and ion loss rates, because charge 
neutrality must be maintained, allows expression of w^ in terms of 
known parameters. Therefore, 


-ek/T nc _ 

£ v w L e ^ hjr r, L 

4ecc 2 e i c 


c _ e 9 JT 

w c = 8 6 


T V 1 * e tp /T 
( -e Jl) p 7 e 

e l V 


(5) 


( 6 ) 


Assuming an argon plasma with T e / T ^ = 1° gives 

* V T . 


w = 0.27 r e 


(7) 




Of course, w c is determined by physics other than that involving the 
interior plasma potential <p^. In fact, if w^ is determined, then 
this equation determines the plasma potential such that the electron 
and ion loss rates are balanced. For example, using w c = r g , then 
gives e cp^ = 1.29 T g , a result which is very close to that experi¬ 
mentally observed in the present work. Thus, cannot be greatly 
different from r g or the predicted plasma potential will not agree 
with the data. Furthermore, in a discussion of the properties of 
the systems as the electrode bias is changed, the electron loss 
along the cusp will be completely dominated by the exponential terms 
in Eq. (l) so the exact value used for is unimportant. Therefore, 
we take w c = r g in the following and write 


( 3 ) 


Now consider the case e cp^/T g » 1. This, of course, results 
from a high electrode bias. In this regime the electrons are 
electrostatically trapped along the cusp, so the only loss process 
for electrons is diffusion across the magnetic field lines. 

Dolan defines a characteristic time for electron loss by 
spatial diffusion to be:^ 
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l£ 


electrode, and is a flux surface that is a distance r frctn the 
’ 0 e 

center of the magnetic cusp. 

Plugging these factors into Eq. ( 9 ) gives 


, 2 2 


/m 

' e 


I « T e V ! 1 B 


(13) 


Let 5 q and B Q be the values of 5^ and B at any arbitrary reference 
point (Fig. 12). By magnetic flux conservation: 


5 


v and T are constant in the main olasma volume because the 
c e 

central plasma is field-free and uniform. We will assume that 
and T g are constant all along By definition dl # = di/ 2 rr r^ and 
B* = B/B q . Then: 


0 -J 


(15) 


X is written:^ 



(16) 






N is the number of magnetic cusps, and P D is the electron gyro-radius 
at the reference point where and B q are evaluated. 

The cross field leak rate for electrons is defined 


where V is the plasma volume. New substituting t gives: 


with S the distance between cusps (« 8 cm) and L, the magnetic cusp 
length. Now is the order of half the electrode gap G. The 
electron gyro-radius at the boundary can be written Pq = 
where fi is electron cyclotron frequency. The exponential is again 
an electrostatic confinement factor with 





Now V/r^ = A the surface area of the plasma: 

A = SL c , ( 19 ) 

with S the distance between cusps (** 8 cm) and L, the magnetic cusp 
length. Now 5 Q is the order of half the electrode gap G. The 
electron gyro-radius at the boundary can be written Pq = v g /ft 
where Q is electron cyclotron frequency. The exponential is again 


an electrostatic confinement factor with 



The cross field leak rate is: 


eLT * « % • t> 


- g e(© - © r )/T 


- g e(© - © )/T 

S k V e W T L c 6 • 


Where in the last equation we define an effective cross-fieId 
leak width, w_, given by 


, can be written, ignoring E field fluctuations, as 


where is the electron-ion collision frequency and is the 

electron-neutral collision frequency, v and \> are estimated 

ei en 


v = U X lO -5 n/T 3 / 2 


t « 5.5 X10 8 p Q T g 5 / 2 


(25) 
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The density n is measured in cm p Q in Torr, in eV, and v 

-1 

in sec 

ft is the cyclotron frequency evaluated midway between the 
cusps. Measuring the field strength along ^ it was found that IB( 
at this point is approximately 20 Gauss. Taking |b| = 0 in the 
interior, we see an average |b| for electron diffusion is 10 Gauss. 
This gives ft = 1.8 X 10^ rad/sec. 

Now taking the ratio of w q and w^, assuming T g ~ 5 eV gives 



This expression with w c /w T normalized to G is graphed in Fig. (l3). 

-U 9 -3 10 -3 

This graph shows that for p Q ~ 10 and 10 cm * n s 10 cm 

then (w c /w t )/G ~ 8. Since 0.6 cm s G s 1.2 cm we see that 

4.8 * w c /w T s 9.6. 

In general, plasma density is determined by a balance of plasma 
production and losses. The creation rate of plasma electrons can be 
written as 01^ where 1^ is the filament discharge current and 6 de¬ 
pends on the cross section for ionizing and nonionizing collisions, 
and on the primary electron confinement time. During an individual 
experiment 01^ was constant, independent of cp £ , therefore 
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PI d 


V eLc + V eLT 


( 27 ^ 


When tp E = 0 and e cp^/T^ ~ 1 means pi d ~ '° e j jC > but when cp v ~ cp p » T e /e 
(cp E ~ cp p will be proved later) then PI^ ~ v eIjT * The last sentence 
refers only to the dominant loss mechanism at a given limit for c? r . 
Equating the two loss rates shows 


n(o) 

4 


v w L e 


e cp p(0 )/T e 


= T~ V w_ L 


(28) 


The density ratio is 



e <P p (0)/T e 


w 




W T 


(29) 


Using the limits of w c /w T means 

1.25 s ^*2.5 . (30) 

These values are in excellent agreement with the ratios observed in 
Fig. (7). 

At sufficiently high neutral pressure or low plasma density 

where electron-neutral collisions dominate ( v e L C / v e i J p) a G / V en aG /P 0 * 

However, if we consider the case of v » v then 
’ ei en 
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^Lc ( 

eLT 


G ( 

V ei 


(31) 


Here the effective cross field leak rate v eLT is proportional to the 

plasma density so the steady state density ratio is no longer 

proportional to the ratio given by Eq. (29). The neutral pressure 

and density at -which the transition from v » v to v » v 

en ei ei en 

depends on the particular plasma species and electron temperature. 
The dependence of the enhancement ratio on the gap width G and the 
discharge current 1^ can he seen as follows. At high electrode bias 
the balance of electron production and loss (when » v gn ) can 
be expressed: 




PI d a °H v e r e L c “ 


(32) 


so the density at high electrode bias goes like a VI^G . At lew- 
electrode bias, when losses along the cusp dominate, the density is 
proportional to 1^ so the ratio of density at high electrode bias to 
zero bias is: 


n 



i. 


(33) 


in agreement with Fig. (7) 











The dependence of the density on electrode bias can be deter¬ 


mined from: 


L [w 


V T , 


■ - v- /T e 


(3^a) 


However, Eq. (jUa) requires a relationship between cp^ and <P E - To 
find this relationship requires equating electron production and loss 


- e <pi/T • g e(® - cp)/T 

-vLfwe p e + we P E e l 

k e c c T J 


■ I % ^vT L = R( y v ■ 


(3Ub) 


where R is a retardation factor that is unknown, and is 

assumed. This equation can be reduced to 


-'VVi -* e (* p -VA e 

e * + p e ^ = 0.2 


(35) 


where W C / W T ~ 8 


has been used. 


Consider first <p_ = 0, but V >0 then Eq. (35) can be written 
t P 


- e <P /T 

e p = 0.27 


(%) 





■where R = 1 because the ions are lost by an unimpeded Bohm flux. 
Therefore 


| (0.27) + J = 0.27 R = 0.3678 , 


or implies R > 1. However, R ^ 1 so cp must increase as some 

P 

function of <P E to keep R ^ 1. Since this function is unknown R 
cannot be calculated, but limits of R can be obtained. 

Assume initially that R = 1 and cp^ > cp £ . Therefore Eq. (39) 
is, solving for <p : 


1 e V T e 


This equation satisfies cp^ > «Pg only if e <P E /T g < 2. Therefore, R 

must be less than 1 when e { P E /T g > 2. Let us now look for R, such 

that = <P and define this R = R . Therefore, Eq. (35) becomes 
b p max 


" e ^/^e 1 

e e + k = 0.27 R 





tPp/T 

°r R may = 3.7(1/8 + e e ). For this experiment 0 * e ? e /T @ 

s 20 (T ~ 3 eV), therefore, 0.46 SR s 1.0. 
e max 

Consider <P^ < «P E , now the only relationsl 
is through R. This regime is the case considered in Eq. (40) 
with cp p replacing <Pg. That is, for R as large as is physically 
possible, cp^ = cp^. We conclude that even if ion losses are retarded 
only to the extent that they can physically balance electron losses, 
we must have cpg ~ <P p . 

The smallest R that is physically reasonable occurs when the 
cusps are plugged by the electrode potential then 


R = / he(tp E - V 7 *! 


( 1 + 2 ) 


T E Y p 


h • 


l E P 

Eq. (35) can be written 


/ 8 V T e + J / e . 0 . w / he( »E - V /T i 

(43) 


Note if g = 1 then h = 0 and vice versa. There cannot be plugging 
for electrons and ions at the same time, Eq. (43) can be solved 
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The solution set of Eq. (43), 
for T g /T^ = 10 and cp^ < cp E , is that eCtPg - cp^) < T g and e(c? E - cp^) is 
largest when e 'Pg/Tg ~ 20 [i.e., e(cp E - cp^) = 0.08 T g ], So the plasma 
potential is slightly less than the electrode potential. We conclude 


equality in the following. 

Let us now consider the density dependence on «P E> Using 
Wj/w c = l/8, and c? E = <p^ » 0, the balance of electron production 
and loss gives 


1 = r v L w [e 

d 4 e c c 




m 


When <P E > 


L c w c <1 *8 ) 


(°)/T 


(45) 


The density ratio is 


- e cp (0)/T 
n 2. e p e 
n(0) " 8 -e <P f /T 


• 1/8 


(46) 


This equation is graphed in Fig. (13). 

Now assume v . » v then Eq. (U4) becomes 
ei en ^ 

eI d“4' , e L c w = (e * ,e/T * * 8nfey* 


(47) 
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where flL « <P » 0. Then: 

E p 

- e q» (0)/T, 

--- (48) 

- e V T e + n 

8n(0) 

This expression is also graphed in Fig. (13)• 

The above considerations show that a theory asstuning the main 
loss mechanism for electrons at large is cross field diffusion 
gives very good qualitative agreement with the experimental results 
of Section III. 

Vacatio n O f . h .vjth q, E the Presence 
of an.Extr.a£tlon_r T pj r d 

Now consider a problem similar to the one just discussed, but 
where one end of e-pot is replaced by a grid or metal plate. The 
calculations discussed here will show good qualitative agreement 
with experimental results presented earlier about simulated ion- 
source operation. 

First, consider the potential of a floating grid. The ion 

current to such a grid with area A g and at a potential <p^ can be 

estimated to be I =enc A/2 provided cp > ^ . I will only have 
S s g P g g 

variation with n, as we showed earlier T is constant, and A is 

e g 

constant. When the grid potential the plasma floating po¬ 

tential, the plasma must supply electron current to balance the ion 
current. At high neutral pressure where the density of the primary 
electrons is much less than the plasma density this balance gives 




. 2S 



1 
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(49) 


or solving for <p - gives 
P * 

T ( - 

cp - cp_ = — fin 2v m i « 5 T /e (50) 

p 1 6 ' m 2TT e 

e 

for an argon plasma. However, at low neutral pressure the primary 

0 

electrons, with energy E^, contribute a current 1 ^ given by: 

I = e n A v cos 9 ( 51 ) 

P P g 


where 


v cos 9 = J <!$■ f(^)v cos 9 


(52) 


The distribution function f(^) is characterized as a shell in 

, .. 16-17 

velocity space 


ftf) 


6 (V - v ) 

_ 2 _ 

Utt v 2 
P 


(53) 


where v is primary electron velocity. Therefore 






Since the right hand side of Eq. (U 9 ) is the same as Eq. ( 57 ) 

requires that - <P g must increase (reducing the plasma electron 

current). This trend in cp - 9 is shown in Fig. (l4). 

P g 

A similar calculation as was used earlier this section can 
employed to calculate the density variation with <P E when the grid 
is in place. With a floating grid the balance of electron produc' 
tion and loss can be written 








This expression ignores primary electrons so is only valid at high 
pressure. This equation also shows that density is reduced by the 
presence of the grid, but the variation with <*> E is not qualitatively 
changed. 

If the grid is grounded, in the same limit as the floating 
grid case, the density varies as: 



( 60 ) 


When <p g < 0 electron current will be less than when tp g = 0. Then 
the last term in the expression for n^ will be less than for the 
last term in Eq. (60). Then the ratio of < 1* If <P g > 0 

the opposite holds and n^/n^ > 1. These arguments are in quali¬ 
tative agreement with the data shown earlier in Fig. (10). 
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V. SECONDARY ELECTRONS 

One consequence of the elevated interior plasma potential is 
the observation of secondary electrons by primary electron bombard¬ 
ment of the chamber walls. These secondary electrons have a 
characteristic signature in the Langmuir probe trace that is clearly 
distinct from the plasma and primary electron contributions. 

The physical picture is a cold electron produced at the wall 
being accelerated through the plasma-wall sheath at the gap and 
forming a shell in velocity space at v = >J 2 e . In ordinary 

multi-dipole devices this velocity is of order v^Te/M e and the 
secondaries merely spread out the Maxwellian plasma. The signature 
left by secondary electrons, when projected onto a planar probe, is 
a straight line starting at cp = 0 (Fig. 15). This straight line is 

-5 

most easily seen when P q < 1 X 10 . At the lower pressure, the 

density of primaries, and therefore secondaries, is a large fraction 
of the total electron density and therefore easily observed. Also, 
at these low pressures the secondaries do not readily thermalize 
with the plasma electrons. 

The number of secondaries produced per incident primary is 
called the yield (6). The 6 wanted is for production at the wall 
and should not include production at the probe. 






The probe current can be written 


I *J vf(v)[i - 6(v )]dv 
V rain 


(61) 


where fi(v^) is the ratio of secondary current to primary current. 

The velocity of electrons at the probe, v w , goes to zero when 
v = v^ n . Energy conservation gives v 2 + 2 e cp^^/m = v^ 2 so 
v^^ = J - 2 e cp^/m where c p^ is the probe bias and e/m is electron 

charge to mass ratio. f(v) is the velocity distribution for a shell 

18 

of monoenergetic electrons projected onto a one dimensional probe 
or 


f(v) « 


1 0 « v * v 

I P 

0 v < v 

' P 


(62) 


where v^ is the primary electron velocity. 6(y^) is known to vary 
with bombarding energy as shown in Fig. l6. For simplicity we assume 


6(v ) » (rp) 


(63) 


where v^ is defined by E u = l/2 Mv^ 2 , and E^ is the secondary 
electron energy such that 8=1 (Fig. 16). Thus, the probe current is 












or 


■{(E 


) (E - e <p ) 


•H 2 ' 


O} (65) 
pr 


where = (l/2) mv^. The primary electron current to the probe 
when the probe is biased at 9^, is the saturation current 



i g 

Hachenberg and Brauer list * 200 V for tantalum, which 

the probe was made from. The assumed form for 6(v^) is a slight 

underestimation, but the maximum error is about 25$ at E = 70 eV. 

P 

A representative probe trace which exhibits all three electron species 
and effects from the probe is shown in Fig. (16). Note the satura¬ 
tion currents of the three components are comparable and that the 
reduction in the apparent primary electron current is significant. 

The fit of the three components is not sensitive to the yield 
coefficient used for secondary electrons emitted at the probe. 

Now that the saturation currents have been corrected they can 
be used to calculate the yield coefficient at the chamber walls. 

Secondaries are created at a rate given by 6n v A where n is the 

P P P P 

density of primaries, v^ is primary electron velocity, and is 

primary loss area. Secondaries are lost either to the walls (A ) or 

s 

by inelastic processes identical to those involving primaries. If 
the secondaries cause an ionization they lose enough energy to be 
trapped and then lost by diffusion (recall Section IV). Balancing 
production and loss gives 








>n v A = n v (A + Vo n ) 
ppp ss s io 


where V is the plasma volume, n Q is the neutral density and o^ is 
the ionization cross section. The probe saturation currents for 
primaries and secondaries, I , are proportional to nv so Eq, (66) 


The term A * r L where r is the primary electron gyroradius, 

P P c P _ 

and L q is total length of the magnetic cusps. Now r^ = . 

Recall from Section II that at the wall B ~ 1 kilogauss so 
r = 3.4 X 10 ~ 5 a/eT\ A typical energy for primaries is ~ 100 eV 
2 

or A = 27 cm . 

P 

The numerator of the bracketed term is not as easily calcu¬ 
lated. It is possible for A g to be as small as r g L q or as large 
as GL, where r g is the secondary electron gyroradius and G is the 
electrode gap, so 6 cm 2 s A g * 800 cm 2 . The ionization term, for 
an argon plasma at neutral pressure p Q , can be evaluated using 
V = 2.2 X 10 1 * cm 5 , o = 2.4 X 10' l6 cm 2 , and n = 3 X 10 l6 p^ (Torr) 


VOn » 1.6 X 10 5 p (Torr) 
i o o 




At p Q = 5 X 10 Torr, this term is .8 cm". Therefore the ioniza¬ 
tion is negligible if A g = 800 cm 2 and non-negligible if A g = 8 cm 2 . 
The extreme estimates of the bracketed term gives 0.33 and 30.0. 

Although E varied for data in this experiment, E is determined by 
p s 

the details of secondary emission and essentially constant so the 

bracketed term can be written a/*/EP . 

Experimental values of I g /l are shown in Fig. (17). The 

values of I /i are comparable to 1.0 so the minimum value of 6 for 
s p 


the general initial trend for 6, where 10 is a representative value 
taken for the numerator of the bracketed term in Eq. (67). 

The above argument has shown that primaries can produce a 
significant number of secondaries. However, the main point was to 
consider the secondary electron contribution to the plasma. The 


electrons degraded in energy by collisions. The ratio of secondary 
production to plasma production is 6 A^/V n Q = 3U 6 when A^ = 27 cm' 
and P Q = 5 X 10 6 Torr. Using the minimum of 0.33 for 6 means that 
6 Ap/Vo^ n Q > 10 so more electrons are being produced by secondary 
emission than by ionization. 

It has been shown that plasma electrons are electrostatically 


confined at the cusps and losses occur by diffusion between the 
cusps. The secondaries are not electrostatically confined unless 
they lose energy. The ratio of density of secondary electrons to 





plasma electrons can be quite high and their probe saturation 
currents can be comparable. Since in conventional multidipole 
devices e w 1, the confinement of plasma and secondary 

electrons is essentially the same and this means neglecting the 
secondary electrons could introduce large errors in calculations 
involving plasma production and loss rates. 

Secondary electrons have not been considered in previous 
multidipole experiments because they are indistinguishable from 
plasma electrons •when wall sheaths are small. In addition, at high 
neutral pressures pr im ary electron densities are relatively low. 

The presence of secondary electrons is not negligible for low neutral 
pressure regimes. Calculations of particle and energy balance in 
multidipole devices that do not take secondaries into account car 
result in significant errors. Secondary electrons are also impor¬ 
tant in interpreting Langmuir probe data. It is apparent that the 
secondary electrons emitted from the probe can lead to substantial 
misjudgement of the primary electron saturation current and conse¬ 
quently the plasma electron temperature (Fig. h ). 
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VI. RELATIVE ION DENSITIES 


This section deals with experiments where the relative density 
of the various ions species are observed. Directly measuring the 
relative ion densities would give clear evidence concerning effec¬ 
tiveness of electrodes in confining the ions. Showing that the 
atomic ion density increases with electrode bias relative to the 
molecular ion density could have a large effect on whether ion 
sources 2 ^ for neutral beam injectors are modified with electrostatic 
plugging. Unfortunately, this measurement could not be performed. 

The relative density measurement was attempted by three methods: 
l) mass spectrometer mounted outside the chamber, 2) particle bursts 
called pseudo-waves, and 3) mass spectrometer mounted inside the 
plasma volume. 


Theory 

The idea behind a mass spectrometer is fairly well known, but 
to refresh the reader's memory the ideas are repeated. Consider a 
particle with mass number A and atomic number Z with energy Ze<p. The 
velocity of this particle is then v = *J2 Ze^/Arn^ where m^ is the 
proton mass. The magnetic field B required to turn the particle in a 
circle of radius r as given by the Lorentz force is A m^ v 2 /r = (Ze/c) 

2 

X vB. There is then a linear relationship between B and e<p. The 
slope is given by (2 nipC 2 /e 2 r 2 ) X (A/Z). In the mass spectrometers 








used here, r was kept constant and ions, assumed monoenergetic, were 
detected after being bent through l80° as B was varied. Flux peaks 
were observed at certain values of B for given cp and A. By plotting 
these critical B values against cp^, the slope of the lines were deter¬ 
mined. Then the ratio of the slopes is equal to the ratios of 
(A/Z) i /(A/Z)j where i and j refer to specific ion species. The slopes 
allow identification of the ion species, but no information about the 
density. The density comes from the area under the peak observed 
with the spectrograph. This area was approximated by using the full 
width at half maximum (FWHM) multiplied by the maximum amplitude of 
the peak. 

The other method involved the use of particle bursts called 

21 

pseudo-waves. Pseudo-waves were first seen by Alexeff et al. in 

1968. These waves were launched as a result of negative potential 

pulses applied to a grid. They explained the launching of pseudo- 

waves from a grid by considering an ion falling through a negative 

potential well, produced by the grid wire. If the pulse was removed 

while the ion was at the bottom of the negative well, then the ion 

will travel on with a velocity v = (2 ecp /Am )‘ where cp is 
max max' p max 

the maximum amplitude of the pulse and a singly charged ion is 
assumed. Likewise, if the particle is somewhere else in the well the 
velocity is v = J 2 ecp/Am^ where cp is the local pulse amplitude. Also, 
positive pulses will work as a pusher of ions away from the grid. 

The reason ion species identification is possible is that 
heavier ions will travel slower so that the ion species will become 
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separate far from the launcher. This characteristic is inherent in 
all waves, however, only pseudo-waves will show a change of velocity 
with pulse amplitude. If you look at the current collected by the 
receiver the time interval between the direct coupled signal and the 
pseudo-wave will decrease as is increased and vice versa. So 

the identification of pseudo-waves is possible and the identification 
of a pseudo-wave for a particular species is possible. The density 
is calculated the same as for the mass spectrometer. 

Experimental Results 

External Mass Spec. 

The external spectrometer was approximately nine inches from 
the plasma anode, and the beam line connecting the spectrometer to 
the chamber was sheathed in n-metal to eliminate the fringing field 
in the tube. The gyroradius was 4 cm, and the maximum magnetic field 
strength was in excess of 2 kilogauss. The electromagnet was 
powered by a programmable power supply, which was swept by a sawtooth 
from an oscilloscope (Fig. l8). 

The spectrographs were recorded on an x-y plotter where the 
x axis was the output of the programmable power supply, and the y 
axis was the output of an electrometer which showed the ion current 
collected. 

A typical spectrograph is shown In Fig. (19) and the magnetic 
field for each ion species is marked on the respective x axis. These 





2 

values for B are used in the B versus V plots to identify the ion 
species (Fig. 20). 

Knowing A/Z, and assuming Z = 1 for a hydrogen plasma, allowe 
identification of H^ + , H 2 + , and H^ + . The quantity of interest is 
the ratio of fluxes for H^ + and H 2 + . The ratio (R) is found by 
measuring the areas under the and H 2 + peaks (Fig. 19). From 
several spectrographs it is possible to see the variation of R with 
<P E (Fig. 20). 

Depending on the choice of p Q and cp^ the graph of R versus 
<P E could be made to do anything. However, taking cp^ - cp^, = 90 V 
and p Q = h X 10 ^ Torr, the curve of the ratio follows very well 
with the curve for plasma electron density (Fig. j). 

There are several difficulties in using an external spectrom¬ 
eter. The beam line from the chamber to the spectrometer places the 
exit port in the middle of a magnetic cusp line. The influence of 
the surface magnetic field could greatly alter the ion extraction 
at low ion velocities and would be different for different masses. 
Also, the long path length allows the possibility of charge exchange 
collisions and a change in the ion spectra. 

Because of these possible errors a method for determining the 
ion species composition needed to be carried out inside the chamber. 

Pseudo-Waves 

The pseudo-waves were launched from a grid with a diameter of 
1.8 cm with 60 lines/inch. The collected current was displayed 
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on an oscilloscope using a high gain plug-in. The pseudo-waves were 
located hy varying the pulse amplitude and locating a wave that 
moved properly. Once the waves were found, photographs were taken 
of the display screen for measurements. 

During the analysis only one pseudo-wave was observed. This 
pseudo-wave had a velocity that an H^ + ion would have if launched 
from the tail of the pulse sheath. A pseudo-wave for a H^ + ion was 
not observed because the velocity of such an ion would be buried in 
the directly coupled signal of the pulse, or if the distance traveled 
by the wave was such that the wave was not buried then damping would 
kill the wave before it could be detected. Due to the pseudo-wave 
not being observed an internal mass spectrometer was built. 

Internal Spectrometer 

The internal spectrometer constructed was based on the design 
22 

of Ehlers et al. This spectrometer could produce a magnetic field 
in excess of 1 kilo-gauss with a fringing field below 100 Gauss. The 
field was produced by two coils with each coil made of 550 turns of 
28 gauge enamel insulated wire. Between the pole faces a copper box 
was used to short out electric fields where the ion current was 
measured. The internal spectrometer had a 1.5 cm gyroradius and a 
two grid entrance aperture. The outer grid was used to keep out 
primary and plasma electrons, the inner grid was used to retard the 
ion energy. The ion energy had to be reduced so the magnetic field, 
lowered by heating, could bend the ion orbit in the required semi¬ 


circle. 



This measurement was also inconclusive due to the poor resolu 
tion of the spectrometer. The poor resolution was a result of the 
fringing field, and that the particle orbit was poorly defined 
because the entrance aperture was poorly collimated. 

The measurement of the fractional increase in could not 
be completed. The methods are theoretically valid but technically 
hard to use. The measurement is important to show that the observed 
increases in the ion density are due to increased confinement. 







VII. SUMMARY 


This thesis has presented experimental results of electro¬ 
static plugging of multidipole cusps. The results show increased 
ion and electron density and increasing plasma potential with 
increasing electrode hias. A model is presented and shown to give 
good qualitative agreement with the experiment. This work has also 
found and measured the contribution of secondary electrons to the 
plasma. 

This work allows other experimentalists the ability to have 
"high" density plasmas at elevated plasma potentials. This freedom 
allows further investigation of secondary electrons, increasing the 
density of a particular ion species over another, or any other 
project requiring this combination of parameters. 
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Figure 3. Typical Langmuir probe characteristic (without secondary electrons). The 
A and B subtract different amounts of primary electron current. 







PROBE BIAS (VOLTS) 


Figure 4 


Semi-log plot of probe current vs. probe voltage (from 
Fig. 3). The effect of incorrectly removing the primary 
electrons is seen to raise the plasma electron tempera- 
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(b) 

Equipotential contours (in volts) determined by an emissive 
probe. 

a. Electrode bias at 2 V, neutral pressure equal to 
8 X 10’5 Torr, n » 2 X 10? cm" 5 . 

b. Electrode bias at 50 V, at 8 X lo" 5 Torr, 
n = 6 X 10° cm" 5 . 
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Figure 9. Variation of experimentally measured loss width with 
electrode bias. Upper figure shows linear relation 
between density and Inverse loss width. 
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Figure 15. Graph of Eq. (26). 
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Figure l6. A representative Langmuir probe trace (labelled DATA). 

The plasma potential was 81.5 V. Data has been displaced 
slightly to allow comparison with the fit to the data. 

The contributions from the plasma electrons l g , pr imar y 

electrons 1^ and secondary electrons I g are separately- 
indicated. The neutral pressure was 5 X 10Torr. 

























MASS SPECTROMETER CIRCUIT 



Figure 19* Schematic of externally mounted mass spectrometer 
experiment. 















Figure 20. Representative spectrograph of a hydrogen plasma. 
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Top three curves show the square of magnetic field for 
maximum flux vs. plasma potential. The bottom graph 
shows variation of density of atomic hydrogen to 
molecular hydrogen with plasma potential when E = 90 
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